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Using positional data from video-microscopy of a two-dimensional colloidal system and from 
simulations of hard discs we determine the wave-vector-dependent normal mode spring constants in 
the supercooled fluid and glassy state, respectively. The emergence of rigidity and the existence of 
a displacement field in amorphous solids is clarified. Continuum elastic theory is used in the limit 
of long wavelengths to analyze the bulk and shear modulus of this amorphous system as a function 
of temperature. The onset of a finite static shear modulus upon cooling marks the fluid/solid 
transition. This provides an opportunity to determine the glass transition temperature Tg in an 
intuitive and precise way. 
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In general, there are different ways to characterize 
supercooled or glassy systems. For instance, a charac- 
teristic change in the thermodynamic properties, e.g. 
volume and enthalpy, upon cooling a fluid may be 
used to define the glass transition temperature Tg [1]. 
Another property often investigated for molecular or 
atomic glasses is the viscosity r]. Given the rapid slowing 
down of the dynamics upon cooling, a system is called a 
glass if its viscosity exceeds lO"'^'' Poise (10"'^^Pa • s) [2]. 
However, this way the transition temperature depends 
on the cooling rate and measurements are often difficult, 
especially in soft matter systems. These systems are 
'softer' up to a factor of 10" (2D) or 10" (3D) 
which interferes with the idea of a universal, viscosity- 
dependent definition of Tg. Since glasses, as crystals, are 
categorized as solids, their mechanical behavior should 
be characterized by a finite shear modulus /i [3]. A 
fiuid, however, lacks such, at least for low frequencies 
[4]. Therefore, we expect to see an abrupt change in 

at the onset of vitrification [3, 5, 6]. Although ^ is 
a macroscopic property, in thermal equilibrium it is 
measurable locally via the equipartition theorem, which 
makes it applicable to soft matter systems [7, 8]. 

In this letter, we analyze the elastic properties for 
both an experimental 2D colloidal glass former and 
for simulation data of binary hard discs. Our anal- 
ysis is based on a quasi- equilibrium description of a 
non-ergodic solid state, which we derive starting from 
fiuid equilibrium \ipon crossing a kinetic glass transi- 
tion. The moduli are derived from thermally excited 
modes in the small wave vector, q ^ 0, limit. A sudden 
rise in the shear modulus marks the onset of vitrification. 

The experimental system consists of a binary mixture 

of super-paramagnetic polystyrene spheres confined to a 
fiat water-air interface [9]. The species A (diameter a = 
4.5 ^m) and B (cr = 2.8 /itm) have a relative concentration 
of ^ = Nb/{Na + Nb) « 45% where Na and Nb are 
the number of particles of both species in the field of 



view. An external magnetic field H, perpendicular to 
the interface, lets us control the particle interactions in 
sif,u. This is expressed by the dimensionless coupling 

parameter 

which acts as an inverse temperature, n denotes the 
2D number density and is computed via Voronoi tes- 
sellation. Xa,b is the susceptibility of species A, B. 
Video microscopy and digital image analysis provides 
the position of individual particles as function of time. 

Our approach is inspired by the one for 2D crys- 
tals [10]. For amorphous solids it requires first an ex- 
perimentally accessible definition of diplacements [6]. 

We start from a fiuid, where a collective mean-squared 
displacement tensor C(q, t) = (Au*(t) Auq(t)) fol- 
lows from time-integrating the velocity field Auq(i) = 
/q dt' Vq(t'). Its Zwanzig-Mori equation of motion 
is determined by stress kernels G(q, f), which reduce 
to the time-dependent rheological moduli for g — >• 
[12]. In a nonergodic, glassy state they take finite 
values at infinite time (index oo), which are related 
by Coo(q) = 2^i^ (Goo(q))"^ Then, the functions 

C(q, t) = \ (Coo(q) — C(q, t)), can be argued to have all 
properties of auto-correlation functions, so that one can 
write them as C(q, t) = (Uq(f)* Uq)^^'''^^ where Uq(t) is 
a displacement fluctuation. The superscript 'glass' indi- 
cates that averaging is done in a restricted phase space set 
by the glassy state. As the so obtained displacement fluc- 
tuations are ergodic, time and ensemble averages agree. 
The former shall be denoted by an overbar so that rl is 
the average position of particle i & [1, A^] during the time 
interval At. Assuming that the particle displacements 
Ui{t) = ri{t) — ¥i are small, the displacement field in an 
amorphous solid follows [12]: 

1 ^ 

Uq(i) = E (2) 



2 



It differs from earlier proposals [13], where coarse- 
grained fields were considered. This completes our 
derivation that the equal-time average of the squared 
displacements is given by elastic coefficients, which are 
frozen-in contributions in g-dependent stress kernels. 
The expression for fl holds as long as At is smaller than 
the relaxation time which was extracted from the 
mean square displacement of the particle coordinates 
for high values of F. Given that in these states the 
mean square displacement does not exceed 1/100 of 
the particles radii within 10** s, we are confident that 
dynamical arrest is experimentally achieved on this 
timescale [14, 15]. 

The Eigenvalues As (the subscript s denotes polariza- 
tion) of the dynamical matrix D^i, are the spring con- 
stants of our system, and can thus be computed from: 



1 



A. 



(3) 



We call the As(q) 'dispersion relation' using quotation 
marks, keeping in mind that the motion is over-damped 
due to the surrounding solvent. By applying this 
method, we implicitly assume only that the system is 
homogeneous and isotropic corresponding to continuum 
elasticity theory. In atomic and molecular glasses, recent 
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FIG. 1. The 'dispersion relation' for different effective temper- 
atures as extracted from the displacements of particles from 
their equilibrium position. Filled and empty symbols rep- 
resent longitudinal and transverse waves, respectively. The 
black vertical line marks the edge of the 'Brillouin Zone' which 
we define as a circle in q space with diameter 27r/a (a being 
the next neighbor distance). 



progress has shown that vibrational modes undergo a 
crossover from propagation to localization due to diffrac- 
tion upon decreasing the wavelength when the local 
structure of the material can not be ignored any more 
[16-18], which may be connected to the 'boson peak' 
(AOP) [19, 20]. In recent colloidal work [21-23], the 
density of states (DOS) was studied and the associated 
modes showed 'swirling motions'. In our appproach, 
DOS and AOP can be studied as functions of frequency 
[24]. 

In Fig. 1 we present the computed 'dispersion relation' 
for different interaction parameters. Filled and empty 
symbols represent spring constants for longitudinal and 
transverse waves, respectively. The growing amplitude 
of the curves for increasing F reflects the expected 
stiffening of the system upon cooling. The black vertical 
line marks the 'Brillouin Zone' edge. Unlike a crystal, 
a disordered (isotropic) system exhibits no symmetry 
points in reciprocal space. Therefore, we define the 
'Brillouin Zone' simply as a circle with diameter 27r/a, 
where a denotes the average inter-particle distance. 
Accordingly, our measured 'dispersion relations' are not 
sensitive to the direction of q, which we confirmed by 
rotating our field of view. In a crystal, the 'dispersion 
relation' scales with the coupling parameter since the 
spring constants are linear in F [10]. Interestingly, we do 
not find this scaling with F in the amorphous solid [12]. 
This hints towards a subtle change in structure since 
particle distances affect the dipolar particle repulsion. 

In high symmetry crystals in 2D and amorphous solids, 
the elasticity tensor C^i/o-t possesses only two indepen- 
dent elements. These two elements can be expressed via 
the Lame coefficients fi and A of continuum elasticity the- 
ory. Following the method described in [7] , we extract the 
Lame coefficients from the 'dispersion relation' as 

a2(2^ + A) 



kr,T 
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lim 
lim 



[g2(|u||(q)|2)si-^] 
[g2(|^,^(q)|2)giassl 



(4) 
(5) 



We chose an intermediate regime (0.8 < qa < 2.0) for 
the extrapolation q — ;> where the data for linear regres- 
sion fits best [7, 11]. In two dimensions, /i gives the shear 
modulus and fi + X the bulk modulus. As we extract the 
elastic constants for q — > 0, we are only accounting for 
long wavelength modes in elastic continuum theory. 

Next, we validate our method by studying finite time 
and size effects, i.e. the frequency and window size de- 
pendence of fJ-{uj). It is known that for high frequen- 
cies, even fluids exhibit a nonzero shear modulus [4]. 
We therefore reduce the time interval At for which av- 
eraged positions are computed. This results in an in- 
creased probing frequency. Fig. 2 clearly indicates the 
expected growth of the shear modulus fi for shorter At for 
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FIG. 2. For high frequencies lu oc 1/At (short trajectories), 
the expected increase in the shear modulus /i (filled circles) 
is observed. On long timescales, the shear modulus does not 
change in the glass (blue dashed line gives best fit) whereas 
it decays for the fluid. For the bulk modulus (empty circles), 
however, we observe a signiflcant decrease upon reaching the 
relaxation time regime. The inset shows a log-log representa- 
tion of the high frequency region where we find a power law 
behavior. 



both liquids (F = 169) and amorphous solids (F = 394). 
These results are consistent with [3]. At the same time, 
we see that for low probing frequencies, the exact value 
of At makes no significant difference for computing the 
shear modulus. This is marked by a plateau of the shear 
modulus in the time domain ranging to values beyond 
the plateau in the mean square displacements where a- 
relaxations are already visible but the intermediate scat- 
tering function did not decay yet to 1/e (At < r^). The 
plateau value of the shear modulus ranges from 10** s up 
to more than 10^ s (which is the longest sampling time). 
Hence, At is set to sa 14400 s for all studied F. 

The bulk modulus of the glassy system shows another 
interesting feature on long timescales (Fig. 2, empty red 
circles). The onset of the observed decay for low frequen- 
cies or long trajectories (lO'^ s) roughly coincides with 
the onset of the relaxation time t„. Structural relaxation 
affects longitudinal modes because their amplitudes are 
extremely small. Transversal modes, however, have 
about 10 times larger amplitudes. Therefore the decay 
of the bulk modulus might indicate structural relaxation 
processes in the a-regime or Mermin- Wagner fluctua- 
tions [25] known from 2D crystals. The inset of Fig. 2 
shows the high frequency behavior of bulk and shear 
modulus. We find a power law dependence for almost 
two decades. The blue lines are guides to the eye with 
exponents of —3/4 and —1/3. 

Furthermore, we examine finite size effects for the 
system at F = 458 in Fig. 3 (a). By varying the 
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FIG. 3. (a) Effect of a reduced field of view on the elastic 
properties. The red dashed lines are guides to the eye. (b) 
Low frequency (long wavelength) mode acquired by diagonal- 
isation of the displacement covariance matrix. The expected 
plane wave character is found. 



window size A/Aq, we reduce the maximum resolvable 
wavelength. A significant change in the extracted bulk 
moduli appears at « 50 % of the largest field of view 
{Aq — 589 /im^). This is illustrated by the red dashed 
line. To a smaller extent, this feature can also be 
found in the shear modulus. To check the nature of 
the long wavelength excitations, the method described 
in [21-23] was used. An eigenvector of the 27V x 2N 



covariance matrix Ci 



=< Ui{t)uj(t) > visualises the 



displacement field of a mode with a frequency given by 
the corresponding eigenvalue. Decomposition of the field 
into the x component for a low frequency mode is shown 
in Fig. 3 (b). The large coherent motion of this proves 
the expected plane wave character. 

The resulting moduli are shown in Fig. 4. As expected, 
the shear modulus /x is zero in the fluid phase (first three 
data points). As we further cool down our sample, we 
find an onset of /i, indicating the beginning of vitrifica- 
tion in which the system becomes rigid (indicated by 
the shading). Its magnitude compares well to a mode 
coupling calculation /i = ISnksT for this system at the 
glass transition [26] . In order to address the question [6] , 
whether the shear modulus jumps to a finite value at the 
glass transition or grows continuously from zero in the 
fluid, we performed Brownian dynamics simulations of a 
binary mixture of hard discs in two dimensions [12, 27]. 
Figure 5 shows data for fj, evaluated according to Eq. (5) 
for different time windows At. The density range spans 
the mode coupling glass transition at (j)c ~ 0.796 [27]. 
The modulus takes finite values as long as the structural 
relaxation time exceeds At. Approaching the glass 
transition from above, /x softens but stays finite, while 
on the fluid side, it vanishes, fi — for At > Tq,. The 
inset of Fig. 5 shows that fi agrees with the long time 
limit of the shear modulus Goo — G{t — >■ oo), supporting 
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FIG. 4. At r « 200, we find an onset of the shear modu- 
lus (red circles). This marks the temperature where the sys- 
tem starts to undergo vitrification. Cooling the system below 
Tg, growing spring constants lead to an increase in stiffness. 
Shading indicates the crossover from fluid (yellow) to glass 
(blue). 



our analysis. 



FIG. 5. Shear moduli jj. from simulations of a glassy mixture 
as function of area fraction (j>. With increasing trajectory 
length (given in units of Dg/a^), fj, in the fluid converges to- 
wards zero. Crossing ^c, /i jumps to a flnite value for At < Tc- 
Inset: Time dependent shear modulus G(t) at = 0.81 calcu- 
lated from the stress-stress autocorrelation function [12]. Its 
plateau value Goo = 124 (magenta horizontal line is given as 
star in the main figure) agrees within statistical error with 
the value oi n = 114 from Eq. (5). 



In conclusion, we have shown that video microscopy 
allows us to study the spatially Fourier-transformed 
displacement field of an amorphous solid. We derived 
its quasi-equilibrium average which, generalizing the 
classical equipartion theorem, is connected to the 
thermal energy and stress correlations. In the long 
wavelength limit we extracted the bulk and shear 
modulus of the system. The discontinuous onset of a 
finite quasi-equilibrium shear modulus separates the 
fluid phase from the amorphous solid. 

P.K. gratefully acknowledge financial support of the 
Young Scholar Fund, University of Konstanz. 
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This supplementary material contains technical details of the theoretical approach (part I) and 
additional simulation results (part II) supporting the conclusions in the main text. Part III contains 
additional experimental information. 



I: DISPLACEMENTS IN GLASS 

In order to define the displacement field in glass, we 
assume the existence of a kinetic glass transition between 
ergodic fluid and nonergodic glass states, and start from 
the ergodic fluid side. Here fluctuations around the equi- 
librium values decorrelate at long times. Crossing the 
glass transition, the states become nonergodic. Fluctu- 
ations (of even time-parity) can not decorrelate and the 
frozen-in correlations describe the glassy state. To pre- 
pare the grounds for capturing local correlations, we re- 
frain from coarse-graining [13]. 



Basic equations 

In a fluid, we consider the velocity autocorrelation 
tensor, K(q, t) = (v*(t) Vq(0)), where for point- 
particles Vq(t) = (l/y^)E»^i e*i-""'(*) hit). Its ini- 
tial value is K(q, 0) = ^^1- Equilibrium averaging is 
done with the Gibbs-Boltzmann distribution. Time de- 
pendence in classical Newtonian systems is given by the 
Liouville operator C, which describes the time deriva- 
tive £ = —idt- In an isotropic system, the tensor of 
velocity correlation functions splits into a longitudinal 
(connected to the compressional) and a transverse ve- 
locity correlator; in two dimensions and with q in x- 
direction: K\\{q,t) = K^^{q,t) = {vx{qx,t)* Vx{qx,0)) 
and K^{q,t) = Kyy{q,t) = {Vy{qx,t)* Vy{qx,0)); off- 
diagonal elements vanish. Similar decompositions hold 
for all following tensors because of isotropy. 

The definition of displacement differences which we 
suggest uses the velocity field: 

AUq(t) = /* dt' Vq(t') . (6) 

Jo 

Because these displacement differences vanish for short 
times and become arbitrarily large in a fluid at long 
times, a well defined quantity is the tensor of collective 
mean-squared displacements: 

C{q,t) = (Au;(t) Auq(t)) = 2 fdt'{t-t')K{q,t'), (7) 

Jo 

where stationarity is used to simplify the two time inte- 
grations to one. 



Equations of motion 

In order to find the displacement correlations at long 
times in a nonergodic state (glass) , the familiar equations 
of motion of the velocity correlation functions are used 
[28]. Newton's equation iCr,{t) = d^r,{t) = ¥.,{i)/m, 
enters, where is the (total) force acting on parti- 
cle i. This leads to the appearance of a memory ker- 
nel M(q,f) = (F* UQ{t) Fq)/{m,kBT) built with the 
force fluctuation Fq = (1///V)Q X)i e^i ""- . The pro- 
jector Q = 1 — V decouples the fluctuations from the 
conserved hydrodynamic ones, density and momentum: 
•p = (5g)q) J-((5g* + Vq) jj^(v*, as the reduced resolvent 

obeys 5t?7g(t) = —i (QCQ) UQ{t). Q also eliminates the 
kinetic part using that velocities are distributed accord- 
ing to a Maxwellian. 

Above equations started from Newtonian dynamics 
and thus apply on times scales much shorter than the 
ones of interest in colloids. Coarse graining to obtain 
equations appropriate for colloids starts from decompos- 
ing the forces F^ on the particles into forces from the sol- 
vent F-°'^ and interparticle forces F?°*. The later arise 
from a potential F?°* = —8^11. We assume that solvent 
and particle forces decouple and that the solvent forces 
are Markovian and can be captured by a friction coef- 
ficient: M(q,i) = ^S{t)l + ^ Gg{t). A prefactor q"^ 
could be pulled in front of the generalized g-dependent 
modulus (a stress tensor autocorrelation function with re- 
duced resolvent) because of Newton's actio-reactio prin- 
ciple J2i Ff"* = 0- Its elements are defined by 

and for vanishing wavevector become the (potential) 
moduli: 

G{t) = lim Gq(i) = ^ ((q • {a ■ q)) , (8b) 

where the microscopic stress tensor element 

and the reduced dynamics simplifies to the full one. Ne- 
glecting the inertial term, the overdamped equations of 
motion of the strain correlators follow: 

C(q,t) + ^^ J' dt' G{q,t-t')C{(i,t') = 2Dotl . 

(9) 
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with the short time diffusion coefficient Dq = ksT/C,. 
To simpUfy the notation, the (wavevector-dependent) 
isothermal compressibihty of the fluid, = Sq/{kBTn), 
which entered via density conservation, is included in the 
moduh with tilde, G(q,f) = G(q,f) + (1/K^)qq. 

Phenomenology of glass 

In an noncrgodic state the time dependent moduli take 
finite values at infinite time: 

G(q,t^oo)^Goo(q), (10) 
which predicts from Eq. (9) 

C(q,t^oo)^Coo(q) = 2^(Goo(q))"' . (11) 

Displacement fluctuations stay below a flnite limit for all 
times. This nonergodic state is a solid one. 

In the glassy state, the relaxation onto the frozen-in 
components can be studied. Experimentally, this re- 
quires the observed state to be deep in the glass in order 
that the structural relaxation is much slower. 

In glass, a displacement fluctuation function C(q, f) 
shall be defined as: 

C(q,t) = ^ (Coo(q)-C(q,t)) . (12) 

Obviously, it satisfies: C(q, i — )■ oo) = 0, and is bounded 
by its initial value: C(q, t = 0) = | Coo(q)- Anticipating 
that one can proof that the C(q, t) are (auto-) correlation 
functions, one thus concludes that they are ergodic. By 
definition, they are stationary. We will assume that all 
familiar relations from ergodicity theory hold [29]. Dis- 
placements in the glass are finite and follow from Eq. 
(6): 

AUq(t) = Uq(t) - Uq = (c^^* - l) Uq (13) 

The C(q, t) can then be written as: 

C(q,t) = (5uq(t)*,5uq)s'--, (14) 

where the superscript 'glass' at the bracket indicates that 
averaging is done in a restricted phase space set by the 
nonergodic glassy state. The displacement fluctuation is 
given by: Juq(tj = Uq(t) - (Uq)«^'^"% where (Uq)s''^« = 

^(|UqP)8lass _ lc^(q). 

Discussion 

The interpretation of the displacement Uq(t) given in 
Eq. (13) remains to be discussed and how to measure 



C(q, t) experimentally in e.g. a colloidal glass. As dis- 
placement fluctuations are ergodic, time and ensemble 
averages agree. Let rl be the temporally averaged posi- 
tion of particle i. From Eq. (13) follows assuming that 
the particle displacements Uj(i) = Vi{t) — rl are small: 

Uq(t) - Uq = / dt'Y, Vi{t') e'*!"-* (1 + 0(q • u,)) 
■^0 i=l 

= f] e^^-^ (u,(t)-u,(0)), (15) 

i=l 

neglecting terms ©(q-u^). Within the first Brillouin zone 
this holds as long as displacements arc miich smaller than 
the average particle separation. Equation (15) leads to 
Eq. (2) in the letter. In this approximation, (uq(/;))s''^^^ = 
0, and the (5-sign for a fiuctuation can be dropped in Eq. 
(14). Prom Eqs. (11,12) follows 

(u;Uq)8'-^ = ^(Goo(q))-' , (16) 

that the equal-time average of the squared displacements 
is given by the inverse of the matrix of elastic constants 
which are the nonergodic contributions in the stress 
memory functions. This central result, used and tested 
in the letter, arises from our definition of displacements 
in glass, which differs from the one presented by Gold- 
hirsch and Goldenberg [13], where \ii{t = Qt) =0 holds 
and the equal-time average vanishes in linear response. 

II: SIMULATION 

We simulated a binary mixture of hard discs under- 
going Brownian motion using the algorithm proposed by 
Scala et al. [31]. The system is made up of = 1000 par- 
ticles, with a diameter ratio of small to big disks ds/db = 
0.7 and equal number concentrations Xs = Xb = l/2 at a 
total packing fraction oi — ^{xgd^^ -\- d^Xb). A de- 
tailed analysis of the structural relaxation close to its 
glass transition can be found in Ref. [27]. 

The dispersion relations and elastic moduli were ob- 
tained as explained in the letter from 10^ equally spaced 
snap-shots along one equilibrated simulation trajectory 
for times up to Af. Following the method of Alder et 
al. [30] we also determined the integrated time depen- 
dent shear modulus 

\ \coHe[0;t] / / 

(17) 

Here the sum runs over all collisions up to time t and 

Ar^j (tc) denotes the relative distance and Apij(tc) the 
momentum transfer of two particles at the collision at 
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FIG. 6. Shear modulus G{t) (symbols) and rjxy(t)/t (solid 
lines) for difTerent packing fractions as depicted in the leg- 
end. Arrows indicate the length of the trajectories, used to 
calculate the values in Fig. 5. 



time tc [32]. The brackets (...) denote the average over 
different simulation runs. The integrated shear modu- 
lus was determined for 600 independent, equilibrated ini- 
tial configurations aX if = 0.81 and for 150 for all other 
packing fractions. Equilibration was assumed when the 
correlation functions became independent on the wait- 
ing time. The differentiation in Eq. (17) was done nu- 
merically, as was the second one to obtain G{t) from 
Vxyit) — Jf^dt'G{t'). Finally, a three point running aver- 
age was performed on the data for G{t). 
Figure 6 shows G{t) and the integrated shear modulus 
Vxyit) divided by t, which is a good indicator of a plateau 
G{t — )■ oo) = Goo. Arrows indicate the lengths of the dif- 
ferent trajectories used in the letter to calculate the shear 
modulus from the average displacement fluctuations; the 
results are shown in Fig. 5 there. Note that the time- 
scale in Fig. 6 is logarithmic: Because of the equidistant 
sampling used to measure (|u(q)p), the final part of the 
time signal dominates. This explains that for (p — 0.80 
the measurement using At = 3472 is below the shorter 
ones. The final a-process already affects this measure- 
ment. 



Ill: EXPERIMENT 

Fig. 7 shows the data of Fig. 1 but normalized with 
the dimensionless interaction strength F. The red 
circles are in the fluid phase while other colors refer 
to the amorphous solid state. They do not collapse 
to a single master curve which hints towards a subtle 



change in structure since particle distances affect the 
dipolar particle repulsion. This change in structure 
does not affect the 2D density n and is barely visible 
in the structure factor. Nevertheless, with increasing 
coupling parameter, particles tend to maximize their 
pair-distances and the wave vector dependent spring 
constant A does not scale linearly with F. 
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FIG. 7. The normalized 'dispersion relations' for different 
effective temperatures in the solid state do not collapse to a 
single curve. The inset shows the extrapolation of the elastic 
constants from the dispersion relation for g — >■ at F = 459. 
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